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Abstract 



Particle emission from primordial black holes (PBHs) has been studied during the past 
three decades, and remains an important testing ground for early- Universe and high-energy 
physics. This thesis gives a review of PBH formation and distribution, the Hawking emission 
process and its later revisions and some of the cosmological consequences associated with 
PBHs. The concept of supersymmetry in quantum mechanics and particle physics is also 
introduced. 

Using the MacGibbon- Webber model, the instantaneous photon (gamma-ray) emission 
spectra from PBHs in the Standard Model and mSUGRA benchmark models were inves- 
tigated using the particle physics event generator Pythia. The point source and diffuse 
gamma-ray spectra for a range of PBH models were also found numerically. 

It was found that the instantaneous photon flux from PBHs in mSUGRA benchmark 
models is generally less than approximately four times greater than in the Standard Model. 
The point source flux from a PBH in its final stages of emission will in mSUGRA benchmark 
models both quantitatively and qualitatively differ from that in the Standard Model. The 
dominant part of the diffuse gamma-ray flux should not be affected by any supersymmetric 
particles. However, for PBHs with initial mass approximately in the range [4 x 10^^ g, 1 x 
10^'' g], the diffuse Standard Model flux at photon energies Eq > 1 GeV will generally be less 
than roughly 50% greater than for mSUGRA benchmark models. The results generally show 
good agreement with previous works, although model differences give some disagreement. 

It is concluded that it is unlikely that the effects found could be detected with the 
upcoming gamma-ray satellite GLAST, due for launch in 2006. 

Key words: primordial black holes, gamma rays, diffuse gamma-ray background, Hawking 
radiation, elementary particle physics, supersymmetry, density perturbations, cosmological 
constraints 
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Chapter 1 

Introduction 



Emission from Primordial Black Holes (PBHs) has been investigated during the past three 
decades, following the works of Hawking [1, 2] which predicted black holes should emit 
quantum particles with a thermal spectrum. 

PBH emission should have a range of effects on cosmological processes, such as entropy 
production, nucleosynthesis and the cosmic microwave background. Photon emission from 
PBHs will contribute to the diffuse gamma-ray background observed. This has been used 
to place constraints on early-Universe conditions. 

The gamma-ray constraints will most certainly improve as the new gamma-ray satellite 
GLAST [3] is launched in 2006. Therefore, it is timely to revise these estimates, and also 
include a possible supersymmetric contribution to the gamma-ray emission spectra. 

MacGibbon and Webber modified the emission picture in 1990-1991 [4, 5, 6] following 
the work by Oliensis and Hill [7] , by assuming only elementary particles are directly emitted 
(rather than composite particles), later fragmenting into hadrons, leptons and gauge bosons. 
This alters the spectrum away from thermality. In this work, we will adopt the standard 
approach of MacGibbon and Webber [4] and apply it to supersymmetric extensions of the 
Standard Model of elementary particle physics. We will find emission spectra for a set of 
mSUGRA benchmark models [8] . 

In Chapter 2, the formation mechanisms and distribution models for PBHs are reviewed. 
In Chapter 3, the emission process is reviewed along with lifetime approximations. Chapter 
4 gives an account of cosmological consequences and constraints, specifically on primordial 
density perturbations, from PBHs. Chapter 5 gives a brief introduction to supersymmctry 
in quantum mechanics and reviews some of the reasons why supersymmctry is interesting 
in particle models. In Chapter 6, expressions for point source and diffuse photon flux 
from PBHs arc derived, the particulars of photon interactions arc presented, an estimate of 
effects from supersymmetric particles on the diffuse flux from PBHs is given and the model 
used presented. The simulation of PBH emission is covered in Chapter 7. The results are 
presented and discussed in Chapter 8, and final conclusions are given in Chapter 9. 

Where convenient, natural units with h = k = c= G= l will be used. 
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Chapter 2 



PBH Formation and 
Distribution 

2.1 Formation Mechanisms 

The original analysis of PBH formation and distribution was made by Carr and Hawking 
[9, 10]. They investigated the formation of PBHs by the gravitational collapse of primordial 
density perturbations [11]. This remains one of the most interesting formation mechanisms, 
especially since it is closely related to very small scales in the early universe. 

However, there arc of course also other possibilities for PBH formation, for example phase 
transitions [12, 13, 14], collapse of cosmic strings [15, 16] and softening of the equation of 
state [17]. 

In our further discussion, we will limit ourselves to PBH formation from the collapse of 
density perturbations. 

2.1.1 Power Spectrum P(k) 

Let (5(x) be the density contrast, i.e. the deviation in energy density relative to the back- 
ground pb, at X, 

pb 

Looking at the Fourier transform of S{x), denoted (5(k), we can write [18] 

(r (k)5(k')) = {2nf6^^\\^ - k>(k) , (2.2) 

where 5^^^ is the 3D Dirac function and P(k) the power spectrum of the density fluctuations. 

This expression assumes the phases are uncorrclatcd, and hence the Fourier transforms are 
uncorrelated. The standard model for P is a power law, i.e. 

P(k) oc jkj" . (2.3) 

Here, n is the so called spectral index. The case n = 1 is the scale-invariant so called 
Harrison- Zel'dovich spectrum. Inflation [18] produces a power spectrum that is close to 
scale-invariant. 



5 



6 



Chapter 2. PBH Formation and Distribution 



2.1.2 Hawking-Carr Model 

Hawking and Carr imposed limits on 5 for black hole formation to be possible. They found, 
using analytical arguments, that a spatial region will form a black hole if the density contrast 
6 in the region fulfills [9, 10] 

1/3<(5<1. (2.4) 

The lower constraint comes from requiring that the size of the region must be greater than 
the Jeans length at the time of collapse, so that gravity can overcome pressure. The upper 

bound comes from requiring that the region not form a disconnected topology, i.e. a separate 
closed universe. The black holes formed in this picture are assumed to have an initial mass 

Mi = j^^^Mh , (2.5) 

where 

^ 1^ (2.6) 

is the horizon mass when the scale of interest crossed the horizon and 7 determines the 
equation of state p = 7p. For radiation domination, we have 7 = 1/3. 



2.1.3 Niemeyer-Jedamzik Model 

A recent development with interesting consequences for PBH formation is the realisation 
that near the threshold of black hole formation, gravitational collapse behaves as a critical 
phenomenon. Numerical simulations by Choptuik [19] and Evans and Coleman [20] suggest 
a scaling relation for black hole formation of the form 

Mi{6) = kMniS - Sc)'' . (2.7) 

Here, fc is a constant dependent on the shape of the fluctuation, 6c the (shape-dependent) 
critical density, 77 the universal shape-independent scaling exponent and Mh the horizon 
mass at PBH formation. 

A valid question to ask is whether this phenomenon will be relevant for black hole 
formation. In general, it will not be relevant for astrophysical black holes, since along with 
other inhibiting effects it would require significant fine tuning of initial conditions. However, 
the effect on PBH formation was further investigated by Nicmcyer and Jedamzik, who found 
that this scaling should be relevant for PBH formation [21]. This is essentially due to the 
fact that for PBH formation, one expects most regions collapsing to black holes to have an 
over-density close to 5c- Further numerical simulations [22] showed that the scaling relation 
holds for PBH formation with nearly-critical initial conditions. They found 

7] « 0.37 (2.8) 

and, significantly, for many fluctuation shapes 

(5c « 0.7, (2.9) 

roughly twice the analytically predicted value of 1/3 (Eq. (2.4)). 
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2.1.4 Epoch of Formation 

Most authors, as will wc, today assume an inflationary scenario, since this solves several 
problems with the standard Big Bang model [18]. Hence any PBHs formed before the 
period of inflation will be diluted away. We consider therefore PBHs formed after inflation, 
typically immediately after the reheating epoch. The time of formation is then related to 
the reheating temperature Tj^h as [23] 



tiRH = 0.3015*-'/'^ , (2.10) 



where g-^, is the degrees of freedom of the constituents in the early Universe (5* ~ 100 in 
the Standard Model and ~ 200 in the MSSM) and mpi the Planck mass. This time is 
easily transformed to the redshift zmH using [18] 

(1/2 \ "2/*^ 

30m Wffo j ^2.11) 

In the Hawking-Carr model, we can also relate the minimum initial mass at a certain 
reheating temperature as [24] 

Mrh ^ l^'/'mp, (^^y\ (2.12) 
where Tpi is the Planck temperature. 



2.2 Spatial Distributions 

The distribution of PBHs is of course closely related to the primordial density perturbations 
(we consider only formation from the collapse of such perturbations). This Section will 
present some models for the PBH mass spectrum, although in the subsequent simulations 
and calculations, these models will not be used (see Section 6.5). 

The present comoving number density of black holes, can be expressed as 

dn dn dMi dn dM, dt dn dt d dM , 

-(1-31717 / -Jirdt). (2.13) 



dM dMi dM dMi dt dM dMi ' dM dt Jq dt 

The two last steps are in practice not very useful, but at least clearly show the dependence 
on dM/dt. In practice, dM/dt most likely yields a separable ODE, from which one easily 
finds Mi as a function of M and t (c.f. Eq. (3.17)). 

2.2.1 Carr Mass Spectrum 

The function dn/dMi, called the initial mass function or sometimes the mass spectrum 
describes the initial number density of the PBHs. The classical result due to Carr is [10] 

which is found using the Hawking-Carr model (see Section 2.1.2), assuming the spectral 
index n = 1. This type of initial mass function is cosmologically irrelevant (see next 
Section). 
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2.2.2 Kim-Lee Meiss Spectrum 

The initial mass function for a general power law power spectrum (in the Hawking-Carr 
model), i.e. P{k) oc fc", was given by Kim-Lee as [25] 

dnKL n + 3 [2 ,,1/2 ,,-5/2 -1 f "f^ \ /oit^\ 

^ = ^\ln^ X exp (^-^ j , (2.15) 

where 7 determines the equation of state (p = jp), pi is the energy density and Mjj^j the 
horizon mass when the PBHs form. This was derived following Carr [10], using the Press- 
Schechter formalism [26] . The mass variance a-j^ , which is roughly the mean square size of 
the density perturbations, will be treated in a following Section. 

The case n = 1 corresponds to a scale-invariant (Harrison-Zel'dovich) spectrum which 
yields a Carr initial mass function, dn/dMi cx ' . As some authors realised, the 
n = 1 spectrum does not yield a significant PBH abundance when normalised to COBE 
observations [27, 28]. Green estimates an initial energy fraction ~ iQ-exw" pgjjg ^g]. 

2.2.3 Niemeyer-Jedamzik Mass Spectrum 

It was found by Green and Liddle [28] that for both power law spectra and flat spectra with 
a spike on a certain scale, in the limit where the number of PBHs formed is small enough to 
satisfy constraints on their abundance, one can assume that all the PBHs form at a single 
horizon mass. 

Using that assumption, one can calculate analytically the initial mass function in the 
Niemeyer-Jedamzik model [21] 



Mi 



driNj ^ Pi / 

dMi ~ V2^rjaHMiMH,i KkMn., 



X 



where fc, 5c and rj are the critical parameters, pi is the c^ncirgy density and MH,i the horizon 
mass when the PBHs form. The mass variance a^j, which is roughly the mean square size 
of the density perturbations, is treated in the following Section. 

2.2.4 Mass Variance 

The horizon mass variance (T(M^f), roughly the mean square size of the density perturba- 
tions, is related to the power spectrum as [30] 

a]j = (t'^{Mh) = / W{kR)P{k)k'^dk , (2.17) 

with R = a~^H~^. W{kR) is a window function picking out scales around k 1/R. 
Using the well-known time evolution of cosmological quantities during radiation and matter 
domination, one finds for power law power spectra [30] 

with (t{Mh,o) = 9.5 x 10"^ from COBE data [30, 31]. The subscript "eq" refers to values 
at matter-radiation equality, and the subscript "0" refers to values at the present epoch. 
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The expression for the mass variance has been reconsidered by Blais, Bringmann, Kiefer 
and Polarski [32, 33], who give a more accurate result, which will not be reproduced here. 
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Chapter 3 

Emission Process 



3.1 Direct Emission 

3.1.1 Hawking Spectrum 

Using semi-classical arguments, Hawking found that black holes should emit quantum me- 
chanical particles [1, 2]. Hawking's expression for the emission spectrum from a black 
hole with angular velocity fl, electric potential $ and surface gravity k, for each degree of 
freedom of particles with energy Q per unit time t is 

(3.1) 



dQdt 2nn (exp (S^^) - (-1)-) ' 



Here, Tg is the absorption probability for the emitted particle, n?i the angular momentum 
and q the particle charge. A black hole is generally characterised by these three basic 
quantities - its mass, charge and angular momentum. This is the essence of the so called 
no hair theorem [34]. For Q = $ = 0, the surface gravity k can be written 



4GM 

If we thus define a temperature T of the black hole. 



(3.2) 



T = = 1.06 C-^] GeV , (3.3) 



SirkGM ■ V M 

we see that when O = $ = the Hawking emission formula mimics a black-body spectrum. 



dQdt 



27rh (exp(§) -(-1)2 



(3.4) 



In the following Sections, we will limit ourselves to non-rotating, uncharged (i.e. Schwarzschild) 
black holes. Page [35] and Gibbons [36] have shown that black holes with spin and/or charge 
will quickly (compared to the times we are interested in) emit particles which carry away 
the spin and charge. Fluctuations will of course still exist, but these have quite a negligible 
effect on the emission rates [4]. 
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Hawking flux 

-2 1 \ 1 1 




2xMQG / ho^ 



Figure 3.1. Hawking flux for s = 0, 1/2, 1 in the relativistic limit. 



A plot of the direct emission spectra for s = 0, 1/2, 1 in the relativistic limit can be seen 
in Figure 3.1. 

The absorption probability Tg can bo further expressed as [4, 37] 

= -2"^ [Q -/i), (3.5) 

where \x is the (rest) mass of the emitted particle and the absorption cross-section. 
The absorption cross-sections have to be numerically computed by solving the Teukolsky 
equation [38, 39]. This was done by Page for s = 1/2,1,2 in 1976-77 [40, 41, 42, 35] (for 
several different scenarios) and for s = by Simkins in 1986 [43]. A plot of the cross 
sections, produced from data in Ref. [4], is shown in Figure 3.2. 



3.1.2 Modified Hawking Spectra 

Hawking's semi-classical derivation has been modified by different authors. We will specif- 
ically look at a recent derivation by Parikh and Wilczek [44]. In contrast to Hawking, the 
authors assume a dynamical geometry, and modeling emission as a tunneling process, arrive 
at 



dQdi 27r?i(exp(87rO(Af-Q/2))- (-1)2-^) • 

We see that the quadratic correction in the exponential becomes important only when Q 
is comparable to M, i.e. in the late stages of evaporation. We then recover the Hawking 
spectrum in the case Q <C Af . These ideas have also recently been used in the context of 
string gravity to arrive at similar corrections [45]. 



Direct Emission 




Figure 3.2. Cross sections for s = 0, 1/2, 1. Fitted using eubic splines from data in [4]. The 
straight line corresponds to the geometric limit blackbody cross section Ug = 277rM^G^/c*. 
The behaviour is similar to an over-damped oscillation around the blackbody limit. For 
parameter values > 0.6 we then approximate with that limit, due to lack of data. The 
accuracy is generally within 5%, increasing with the parameter value. However, as can be 
seen, for parameter values slightly greater than 0.6 the accuracy for s = 1 is only within 
~ 10%. 
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3.2 Fragmentation and Hadronisation 

In the model proposed by Oliensis and Hill [7] and expanded upon by MacGibbon and 
Webber [4, 5], all fundamental particles with rest mass less than or equal to the black hole 
temperature are emitted (of course, occasionally heavier particles are emitted, but as seen in 
Hawking's expression they are exponentially damped). When the black hole temperature 
reaches the QCD confinement scale Kqcd = 250 — 300 MeV, quark and gluon jets are 
emitted. These fragment and eventually produce hadrons and other particles. 

The authors argue that almost all emitted particles are much smaller than the black 
hole, since the de Broglie wavelength of most emitted particles satisfies 



AdB < 0.04 X 10"^^ I — — - I cm <C rbh 



T \" 2GM 



IGeVy c- 



2 



Furthermore, successive emissions should not interact, since for most particles 

1 20 

A.. = ^«At«-, (3.8) 

where Q is the particle energy and is the time between successive emissions. Oliensis 
finds numerically that 99.9% of the emitted particles should satisfy this. 

Likewise, the previously emitted particles should not affect the emission process itself, 
since 

r6/.~^«At. (3.9) 

There should also be no strong interactions between successively emitted particles. 
Lastly, gravitational effects on lifetimes, etc. should not have a significant effect on the 
final spectrum. 

In summary, no short or long range forces should significantly affect the emission process 
before fragmentation and hadronisation, which means the emission is similar to that of 
e"'"e~-annihilation events. 

However, it should be noted that these ideas have been challenged by Heckler. He 
proposed in 1997 that QED and QCD photospheres (or perhaps better termed photosphere 
and chromosphere as for the Sun) may form around the black hole above a certain black 
hole temperature when bremsstrahlung and pair-production becomes non-negligible [46, 
47]. These ideas remain controversial, but seem at least partly supported by numerical 
simulations by Cline, Mostoslavsky and Servant [48] and Daghigh and Kapusta [49]. In 
light of the uncertainties associated with photosphere efii'ects, we will neglect photosphere 
effects in this work. Barrau et al. [37] use a parameterisation based on Ref. [48] to account 
for photosphere effects, but this method seems dubious. 

To find the actual emitted spectra after fragmentation, one has to convolve the direct 
(i.e. Hawking) emission spectrum with a function describing what final particles are created 
by the fragmentation of an initial particle. For a particle species X we can write 

dEdt ^Jn-^ 'o^%f^^^fQ \ _t_,\2s\ dE ^' ^ ^ 



Q=E 2Trh (exp - (-1)2-) 



where we sum over all particle species j and weight with their degrees of freedom aj. The 
fragmentation function dgjx{Q,E)/dE is the number of X particles with energy between 
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E and E + dE created by the particle (jet) j with energy Q. Here we can (and should) 
include any decay products, so for most simulation purposes one may readily decay all but 
particles stable on astrophysical time scales, i.e. 7, pp, e^, vV. This list should also contain 
the lightest supersymmetric particle (LSP) in the case of a supersymmetric model (with 
R-parity conserved). 

One should specifically note that since we have for pions m,r ~ 135—140 MeV < Aqcd, 
in the (approximate) energy interval [m„,AQCD] the black hole will directly emit pions. 



3.3 Mass Loss Rate and Lifetime 

From the emission rate we can find the black hole mass loss rate as [5] 

where we sum the direct emission rate, Eq. (3.4), over all particle species ) with rest mass 
/Xj and aj degrees of freedom. Using the mass loss rate, the lifetime r of a black hole is 



-j 



'"^^ dM 

i^y'dM, (3.12) 



where the inverse mass loss rate is integrated from the initial black hole mass Mi to the 
Planck mass mpi. The latter limit is an approximation to deal with the fact that our 
derivations, particularly Hawking's, break down close to the Planck scale. Uncertainties in 
the final mass will have a negligible effect on the lifetime (see e.g. Eq. (3.17)). 

Similarly, we find the initial mass Mj of a black hole formed at f = 0, that would just 
have evaporated at time tf to be 

M. = 0dt. (3.13) 

If we assume all PBHs are formed very close to the Big Bang, we can replace tf by t„, the 
age of the universe, to get the initial mass of PBHs that would just have evaporated 
today. The generally used value is M* = 5 x 10^** g [50]. 
Far from mass thresholds, one can assume that [50, 5] 

^«-$(M)M-2, (3.14) 

where $(M) then is a function of the number of emitted particle species. This function 
generally has to be calculated numerically, and the values reported by MacGibbon and 
Webber [4, 5] are 

^ = -5.34xlO^V(M)(^)'' gs-i, (3.15) 

where /(M) in the relativistic regimes, i.e. far from mass thresholds, is 

/(M) = 0.267rfs=o + 0.060ds=i + 0.007rfs=2 + 0.020^^=3/2 + 

+0.147Cv2 + 0.142C?/2 . (3-16) 

where the factors d count the degrees of freedom for particles with spin and charge according 
to the sub- and superscripts, that are being emitted by the black hole (this is how the M- 
dependence comes in). The s = 3/2 contribution is for hypothetical gravitino emission. The 
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Mass loss rate function f(M) 
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Figure 3.3. Mass loss rate coefBcient function /(M). 



error in counting quarks as having charge ±e or charged elementary particles as uncharged 
is < 5% [4]. This function is plotted in Figure 3.3 for the Standard Model and mSUGRA 
Benchmark B. MacGibbon also made approximations to account for the effect of mass 
thresholds [5]. 

Using the assumptions (3.14)-(3.16), Eqs. (3.11) and (3.12) can be solved for M and r 
respectively, to find 

M ={Mf- 3$i)^/^ (3.17) 

and 

when $(M) is roughly constant. MacGibbon estimates the black hole lifetime to be [5] 

T « 6.24 X 10^27 (^Ml^ f{M,)-^ s, (3.19) 

with /(M) given by Eq. (3.15). This is naturally a slight overestimate of the lifetime, since 
the increase in degrees of freedom as the temperature rises is neglected. However, since the 
black hole spends most of its time near Mi the relative error is small. This estimate of the 
PBH lifetime is shown in Figure 3.4, where the effect of the aforementioned error can be 
seen in that the mSUGRA Benchmark B lifetimes start to differ from the Standard Model 
only at masses lower than approximately 10^^ g. 

These expressions for the mass and lifetime hint at why PBHs are particularly interest- 
ing, i.e. why we do not consider emission from black holes formed from the collapse of stars 
at later epochs. A black hole with M = Mq « 2 x 10^"^ g would have an initial temperature 
Ti ^ 10^^^ GeV and a lifetime that is ~ 10^^ times the age of the universe. For all practical 
purposes, these black holes will thus emit nothing. 
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Figure 3.4. Black hole lifetime as a function of mass 
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Chapter 4 

Cosmological Consequences 



4.1 Effects on Cosmological Phenomena 

Hawking emission provides a mechanism for PBHs to affect the early Universe. The emitted 
particles can interact with radiation and matter in the Universe and thus affect many of the 
early- Universe phenomena: the CMBR [51, 52], entropy production [52], baryogenesis [53] 
and nucleosynthesis [54, 55, 56, 57] to name some. Recently is was also realiscid that PBHs 
are potentially a great probe of the early Universe with a varying gravitational constant 
[58]. 

4.2 Physics Constraints 

PBHs can be used to probe the primordial mass distribution, hence the spectral index n 
in a power law power spectrum. Page and Hawking required that the emitted gamma rays 
from PBHs should not exceed the observed diffuse gamma-ray background (DGB), and thus 
found a constraint on the abundance of PBHs [50] . This is known as the Page-Hawking 
bound. This constraint has later been updated by MacGibbon and Carr [59, 60] using 
data from EGRET [61] and COMPTEL [62]. Kim, Lee and MacGibbon [24] also use the 
normalisation to the DGB to constrain the primordial density fluctuations in the case of a 
Kim-Lee initial mass function. They find n < 1.23 — 1.25 for the spectral index. 

Another popular approach is to require that some critical energy fraction Q < 1, typically 
the total energy in PBHs Q.pbh- This kind of approach was used by Carr, Gilbert and 
Lidsey to find constraints on the spectral index n [63]. They find n < 1.3 — 1.4. Green 
and Liddle have corrected errors in their results and find n < 1.25 — 1.3 [30]. Kim and 
Lee performed an analysis using the Kim-Lee mass spectrum and obtained n < 1.43 — 1.48 
[25]. Kribs, Leibovich and Rothstein used the Niemeyer-Jedamzik mass spectrum to obtain 
n < 1.2 — 1.45 [64] (they also use the DGB bound). All these calculations have been done 
using normalisation to COBE anisotropy measurements [31, 65]. 

Other methods of constraining n are also used in the aforementioned and other articles. 
Reviews of these and other cosmological constraints from PBHs are given in Refs. [66] and 
[60]. 

Currently, the best constraints on n come from CMBR measurements such as COBE 

[31, 65] and WMAP [67]. However, it is important to note that n may vary with fc, and 
PBH constraints probe much smaller scales than CMBR measurements do. Consequently, 
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the PBH constraints provide information that CMBR measurements do not, especially in 
the case that P{k) is not a power law but a spike or step [33]. 



Chapter 5 



Supersymmetry 



5.1 Supersymmetry in Quantum Mechanics 

To understand the basic concept behind supersymmetry (SUSY), let us start from the 
quantum mechanical Hamiltonian for a single particle: 

We now assume that the potential Vi (x) is unknown, but we know the ground state wave 
function tpoix) which we will assign the energy eigenvalue (no loss of generality) . Inserting 
the ground state wave function into the Schrodinger equation for the above Hamiltonian we 
find 

V.{x) = f^. (5.2) 
^ ' 2mVo(a;) ^ ' 

If we now set out to write the Hamiltonian Hi in a standard form, i. e. Hi = A^A, we can 
write 

A = -^4- + W{x), At = --L-^^ + W{x) , (5.3) 

where W{x) is the so called superpotential. Looking at the original expression, we can now 
write 

Vi{x) = W^{x) - -^W'{x) , (5.4) 

V2to 

which is called the Riccati equation. Using the fact that Atpo iJi^o = we find the 
solution for W in terms of the ground state wave function 

n (5.5) 
Let us now define the operator H2 = AA'^ . We then find 
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V2{x) = W^ix) + ^.w'{x). (5.6) 

We call Vi {x) and V2 {x) supersymmetric partner potentials. Using the Schrodinger equation 
for the two Hamiltonians, one finds that the eigenvalues and eigenfunctions are related as 

E^:^=Ei%Ei'^=0, (5.7) 

= (E'^^lr^'M^ , (5.8) 

i^i% = {E(^^)-'/'A^^:K (5.9) 

Both sets of eigenvalues are positive semi-definite, and the energy levels are degenerate, 
except for Hi which has an "extra" level at zero energy. We also see that one effect of 
the operators A and A'^ is to "jump between the potentials". One can now define a SUSY 
matrix Hamiltonian 



and the operators 



\ ^+ / At 



A - )• (5-11) 



These three operators now form the closed superalgebra sZ(l/l) with commutation and 
anticommutation relations 

[H,Q] = [H,Q^] = 
{Q,Qt} = H 

{Q,Q} = {Q\Q^} = (5.12) 

The so called supercharges Q, commute with H, hence the degeneracy in the spectra of 
Hi and H2. The supercharges can be interpreted as operators that exchange bosonic and 
fermionic degrees of freedom. 

In broader terms, what we have seen is that given a particle (ground state wave function), 
we may construct another particle such that they are both described by the same superpo- 
tential, and bosonic degrees of freedom in one particle correspond to fermionic degrees of 
freedom in the other. In a sense they are both "components" of a "superparticle" . 



5.2 Supersymmetry and Particle Physics 
5.2.1 Background 

The ideas laid out in the preceding Section form the basis of supersymmetric extensions of 
the Standard Model. The main motivations for considering supersymmetry are: 

• solution of hierarchy and naturalness problems 

• connection between gravity and other fundamental interactions 

• can learn about high-energy physics in general 

• rich experimental consequences 
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The hierarchy and naturalness problems are intimately related. Suppose we wish to write 
down a Lagrangian which accounts correctly for physics up to some energy scale A much 
larger than the electro-weak scale E^gak ~ 200 GeV. Choosing to break electro-weak sym- 
metry by the Higgs mechanism, one finds that the masses of scalar excitations of the 
field will be of order A. This is due to radiative corrections, the mass corrections being 
^''^Icaiar ~ 9^^- Specifically, the Higgs mass would become huge, destroying the hierarchy 
between the electroweak and A scales, giving us the hierarchy problem. The simplest way to 
avoid this problem is to fine-tune the coupling constants to an accuracy of E'^eakl^^ ~ 
(in the case A = Equt ~ 10^^ GeV). This is the naturalness problem. Supersymmetry re- 
moves or significantly improves on these problems by introducing new corrections from the 
new particles, which reduces the fine-tuning requirement to at most \og{h./ Eyj^ak) and com- 
pletely removes it in some cases. Since we have not observed supersymmetric particles with 
the same masses as their Standard Model partners, supersymmetry can not be an exact 
symmetry. We thus need a mechanism for spontaneously breaking supersymmetry. The 
symmetry breaking imposes approximate limits on the mass splitting between a particle 
and its supersymmetric partner: if the Higgs boson mass correction 5m^^ ^ — '77,Ji')j 

where tob is the boson mass and iriF the fermion mass of a particle-supersymmetric partner 
pair, is to be of the order of the electroweak gauge-boson masses or less, the mass splittings 
should not be much greater than the electroweak scale. 

On a side note, there are other models for dealing with the naturalness problem, for in- 
stance the Technicolor approach where fundamental scalars are composites of new fermions. 
This approach is however technically difficult. 

Local supersymmetry, i.e. where transformations are allowed to depend on position 
and time, and gravity are closely related. Loosely speaking, local supersymmetry gives 
local Poincarc symmetry, which is the basis for general relativity. Because of this relation, 
local supersymmetry is called supergravity (SUGRA). Initially, it was thought that this 
supergravity could hold as a successful theory of everything. However, theoretical consider- 
ations have shown that it most likely would only be an effective theory up to some energy 
scale. Nevertheless, it remains important since it turns out to be the proper framework for 
spontaneous breaking of supersymmetry (which is required) . 

There are rich experimental consequences from supersymmetric particle models, so it 
is possible to test the ideas. One important prediction is the existence of stable particles, 
which could explain at least part of the cold dark matter problem in cosmology. 

Even if supersymmetry in the end turns out not to be a correct description of Nature, we 
can learn a lot about how high-energy physics should behave from studying these theories. 

There are also aesthetic and philosophical arguments for considering supersymmetry. 
For instance, Nature has shown it "likes" gauge symmetries, and supersymmetry would be 
the next natural such symmetry. 

A review of supersymmetric matter and observational aspects can be found in ref . [68] . 
Ref. [69] gives a general introduction to supersymmetry in quantum mechanics. A com- 
prehensive introduction to supersymmetry and supergravity is given in Ref. [70]. The 
discussion above is based on those references. 

5.2.2 Minimal Supersymmetric Standard Model and mSUGRA 
Definition 

The Minimal Supersymmetric Standard Model (MSSM) is, as the name implies, the stan- 
dard supersymmetric model. It contains all the known particles in the Standard Model plus 
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an extra Higgs isospin doublet (that is theoretically required) together with their supersym- 
metric partners. The interactions of the model are those allowed by SU{3) x SU{2) x [/(I) 
symmetry and renormalisability. The MSSM, as do most supersymmetric models, requires 
R-parity, defined as 



to be multiplicatively conserved {B is baryon number, L lepton number and 5 spin). Stan- 
dard Model particles have R = 1 and supersymmetric particles have i? = — 1. This places 
constraints on the allowed processes. For instance, we see that there must be a lightest 
supersymmetric particle that is stable (the LSP), since it has no allowed decay mode. Also, 
supersymmetric particles produced from Standard Model particles will be pair-produced. 

The most general MSSM contains 63 parameters, disregarding the usual Standard Model 
parameters. For practical purposes, some form of restrictions are usually imposed, limiting 
the number of free parameters. 

We will use a parameterisation known as mSUGRA (also UHM, denoting "Unification 
of Higgs Masses" or CMSSM, "Constrained MSSM") which consists of five parameters, the 
common scalar mass mo, the common gaugino mass mi/2, the sign of the higgsino mass 
parameter sgn(/i), the ratio of Higgs vacuum expectation values tan/3 and the common 
trilinear coupling A (see next Section for an introduction to these particles). This model 
assumes that the supersymmetry breaking is gravity-mediated, and in this framework the 
masses toq and TOi/2 are assumed universal at some GUT scale. This parameter space 
has been investigated by Ellis et al. [8] applying various experimental and theoretical 
constraints. They give a set of benchmark models (parameter sets) that represent parts of 
the parameter space that are particularly interesting in light of the constraints. We will use 
these benchmark models in this work, referring to them as "mSUGRA Benchmark A-M" 
(see also Appendix A). 



Particle Content 



For each ordinary Standard Model particle we have a supersymmetric partner, or sparticle. 
The naming convention is such that fermions gain a prefix 5 and bosons a postfix ino. For 
instance, quarks have squark partners and photons photino partners. Table 5.1 should make 
the notation conventions clear. 



R= (_l)3(B-i)+2S' 



(5.13) 
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Particles 


Sparticles 


Name 


Symbol 


Spin 


Name 


Symbol 


Spin 


Charged 


£^ 


1/2 


Charged 







leptons 






slcptons 






Neutrinos ^ 




1/2 


Sncutrinos 







Up quarks 




1/2 


Up squarks 


~„ ~u 

Il.b- Il.r 





Down (jUiiiks 




1/2 


Down stjUfnlvS 


iL.R IL.R. 





Gluon 


g 


1 


Gluino 


g 


1/2 

/ 


W-bosons 




1 








Charged Higgs 







Charginos 




1/2 


boson 












Photon 


7 


1 








Z-boson 




1 








Light scalar 


h? 











Higgs boson 






Neutralinos 


Al,2,3,4 


1/2 


Heavy scalar 


H° 









Higgs boson 












Pseudoscalar 


AO 











Higgs boson 












Graviton 


G 


2 


Gravitino 


G 


3/2 



* Since neutrinos actually have a small mass, we should expect there to 
be a right-handed version also, although this has never been observed. 
Consequently, the sparticle side should in principle also contain corre- 
sponding degrees of freedom. 



Table 5.1. Standard Model particles and their SUSY partner sparticles. 

Table 5.1 lists charginos and neutralinos, which are the mass eigenstates one finds by 

diagonalising the Hamiltonian. Hence these states are the natural ones to use in most 
physical contexts. They can be written as superpositions of photinos, zinos, winos and 
higgsinos, i.e. 

Xf = Mf,W^+Mt,H^, (5.14) 
X° = 7Vii7 + A^i2^° + A^i3^? + A^i4^2°> (5.15) 

where the Ms and Ns are the usual expansion coefficients. 

In this work, we use two different meanings when referring to "the Standard Model". 
When considering "the Standard Model" in its own right, e.g. in contrast to a supersym- 
metric particle model, "the Standard Model" will be understood not to contain any Higgs 
particles. When considering a supersymmetric particle model in its own right, Higgs bosons 
will be looked upon as "Standard Model" particles, although formally not the same "Stan- 
dard Model" as mentioned in the previous sentence (one should perhaps rather use "the 
Extended Standard Model" in this case). The notation of Table 5.1 might confuse this 
issue, so the reader is advised to take special note of these conventions. 
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Chapter 6 

Photon Spectra from PBHs 



6.1 Point Source Flux 

Choosing our position in the universe as the origin of a spherical coordinate system (r, 0,(j)), 
we can express the (unattenuated) photon flux at Earth from a single PBH at r as 

dNf{r,M,E) _ 1 (PN^{M,E) 

dEdtdA 47rr2 dEdt ' ^ ' ^ 

The symbol ® is the astronomical symbol for Earth. In principle, we should of course include 

redshift dependence and effects from interactions. However, since we can only expect to 
detect evaporating PBHs within ~ 1 pc [29] we neglect those effects. 



6.2 Diffuse Flux 

The diffuse gamma-ray spectrum reaching Earth will in general be dependent on dM/dt, 
dn/dMi, the cosmological model used and the interactions of the gamma rays while traveling 
towards Earth. These are suitably described using the optical depth T{E,t), which is the 
effective number of interactions. 

If we consider an isotropic distribution dn/dM of PBHs, one can easily write down the 
flux of photons with energy between E and E + dE within the solid angle fio from PBHs 
with mass between M and M + dM as 



dN®{M,E,no ) _ C'"^-^'' dNf{r,M,E)dn{M) 
dEdtdAdM 



r dNf{r,M,E)dn{M) 

Jo Jo dEdtdA dM ^^■^> 



Since we assume an isotropic PBH distribution we can perform the solid angle integral to 
obtain 

dN^iM, E, no) _ no d-'N^iM, E) dn{M)^^ 



dEdtdAdM 4tt Jo dEdt dM 

where we have made use of the expression for the flux from a single PBH, Eq. (6.1). To 
turn this into a useful expression, we note that since we are dealing with photons, we have 
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dr = —cdt (the minus sign since the photons are travehng inwards in the chosen coordinate 
system). We can thus write 

rfjV,?(M,g) _ c /-^o ffj^j,.dn{M) 
dEdtdMMdn, -^l . ^-^(^'^)^Ar^*^ ' ^^-^^ 

where for simplicity we introduce 

MM,E)= ^^^^ \e=e', (6.5) 

the instantaneous photon flux at energy E' from a PBH with mass M. We should also 
integrate over M to get the total spectrum, 

dN®{E) c /•'" /■^'""^^*'^ , j.,dn{Mi) 



dEdtdAdflo 4?: J^^ 



/ / f-y{Me,ap,E)—^dMidtu (6.6) 

Jt^,„ JM,{ti) d-^^i 



where we have changed variables to the initial PBH mass Mj. The limits are given by the 
range of masses that are emitting at time t\, which will be [M^,{t\),M^ax{t\)\ {J^max{ii) 
is the heaviest existing PBH at time ti). The quantity M^yapiMi,'!'') is the mass at time t 
of a PBH with initial mass Mj, evaluated according to Eq. (3.17). This expression neglects 
effects from the expansion of the universe and degradation due to interactions (the optical 
depth). Let a(t) be the FLRW scale factor (denoting aj = a{tj)). The relevant expansion 
effects are photon redshift and the decrease in volume density. Including these and the 
optical depth, we thus get 



dJ_ ^ dN^jEp) ^ c_ r^'^^^i^^) ao /oi 

dEo ~ dEodtdflodA Jtrr^ir. J M.{ti) «i K^-i 



^-t{Eq,zi) ^ 



X A(Me„ap, (1 + zi)Eo)^^^:^^dMidti , (6.7) 

where E = (l + zi)Eq was used to rewrite in terms of the observed photon energy Eq and 
Zi = z{ti). In keeping with convention, subscripts i denote the time of PBH formation. We 
may also rewrite this expression completely in terms of the redshift z using the definition of 
redshift, 1 + Zj = ao/aj, and the following expression (for derivation, see for instance Ref. 
[18]). 

^ = -Ho\l + z)-'[{l + zf{l + Qmz) - z{2 + 0)Oa + fifl(l + z)^'/^ . (6.8) 
dz 

Here we clearly see the dependence on the cosmological parameters. The radiation contribu- 
tion flji is only relevant for very high redshifts in the radiation-dominated epoch (z > 1100). 

This derivation of the diffuse flux follows those of Hawking and Page [50] and Kim, Lee 
and MacGibbon [24]. 



6.3 Optical Depth 

The optical depth t{Eo,z) has been studied by Zdziarski and Svensson [71]. We will 
look at the processes they consider in the following subsections. Since we are considering a 
supersymmetric scenario, one also has to consider potential contributions from such particles 
to the optical depth. However, any such sparticles would quickly decay to ordinary particles 
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and LSPs. The LSP-photon cross sections are very small, so we should not expect any 
difference in the optical depth due to sparticles. 

For an introduction to the notion of optical depth, see e.g. Ref. [72] . 

In the following, we will use the same parameterisation as Zdziarski and Svensson, i.e. 

T..r ^ (6.10) 

ho ^ ^, (6.11) 

Oo.i = ^, (6.12) 
Go = 4.55 X 10-^°T2.7 . (6.13) 
Furthermore, the formulae given below assume 0,tot = 1- 

6.3.1 Photoionisation 

Neutral H and He atoms may become photoionised by photons with energies larger than 
the ionisation energies, 13.6 eV and 24.6 oV respectively (7A A+e~). 

Below the reionisation redshift Zri, the optical depth chic to photoionisation is essentially 
zero. With that assumption, the optical depth due to photoionisation, Tpi, takes the form 

Tpi = 1.1 X 10-i°Oo.i/i5oeo [(1 + ^h)"'-^ - (1 + z)-^-^] . (6.14) 
The expression for Tpi is a numerical fit valid for photon energies E > 250 eV. 

6.3.2 Compton Scattering 

We divide the effects of Compton scattering (76 je) in a scattering loss and an energy 
loss component. 

Scattering Loss 

The optical depth Tcs due to the scattering component is given exactly by 

Tee = 2.27 X 10-3i7o.i/i5oeo {Fes [eo (1 + z)] - Fes (eo)} , (6.15) 

where 

_ (l + e)^-l/3 3e + 7/4 4 

Pcsie) - ^372 ln(l + 2e)-^^-^e / 

llarctanr(2e)i/21 

— -■ 6.16 

3 X 25/2 ' 

Energy Loss 

The optical depth Tee due to the energy loss component is given exactly by 

Tee = 1-14 X 10-3f2o.l/l50eo ''^^ {^ce [^0 (1 + z)] - Fee (cq)} , (6.17) 



where 



^ , , 2(1 + e)2 , ^ ^ 272g-'' + 470e2 + 261e + 48 
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+ 



Tarctan [(2e)V2] 
2^72 



(6.18) 



6.3.3 Photon-Matter Pair Production 



Pair production from photon-matter interaction is primarily in the form of photon-atom 
pair production (7A jAe~^e^) for z < Zrec and photon-ionised matter pair production 
(76 ee~^,^Z — > Ze+e~) for z > Zrec, where z^ec ~ 1300 is the redshift of recombination. 



Atoms 

The optical depth Tapp due to photon-atom pair production is given by 

Tapp = 9.3 X 10-6fio.l/l50eo "^^^ {Fapp [eo (1 + Z)] - Fapp (eo)} 



(6.19) 



where 



F.pp{e) = e3/^ln(^)- 1650eV2 + 



,825 -Fe 
+2 X 825^/^ arctan 



V825/ 



(6.20) 



The expression for Tapp is found using fitted numerical data, and is valid for (1 -|- z)eo > 6. 
Ionised Matter 

The optical depth Tipp due to photon-ionised matter pair production is given by 



'ipp 



where 



1.15 X lO-^flo.ihno {(1 + zf^ Fipp [eo (1 + z)] - F^pp (eo)} 

109 2 



Fipp{e) = ln(2e) - 
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(6.21) 
(6.22) 



This expression for Tipp is valid for (1 -|- z)eo 1. 



(6.23) 



6.3.4 Photon-Photon Scattering 

Photon-photon scattering (77 — > 77) gives an optical depth r^g 

r^, = 2.44 X IQ-^^h^^Tl^el [(1 + zf^'^ - 
This expression is valid for (1 -|- 2:)^eoOo < 0.1. 

6.3.5 Photon-Photon Pair Production 

Single Photon Pair Production 

The optical depth T^spp due to single photon pair production (77 e^e~) is divided in 
two limits. For (1 -|- z)'^eoQo <C 1 we have 



3.83 X IQ^h^^Tlj (TreoOo) 
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where 



X {(1 + zf^ F,ppi [eo (l + z), Oo (1 + z)] - F,„i (eo, @o)} , (6.24) 



Fs„i (e, 6) = A + e- , (6.25) 



and for (1 + z)-^eoQo 1 we have 

4^2 



Tyspp — 3.83 X 10 



where 



3eo0o 

{Pspph (eo, eo) - (1 + z)"'/' F,pph [eo (1 + 2) , } , (6.26) 
J^«ppft(e,e) = ln(25.5ee) . (6.27) 



Double Photon Pair Production 

The optical depth r^dpp due to double photon pair production (77 — > e+e~e+e~) is given 
by 



■^dpp = 47.5h-„'Tl, (l + z)3/2_i . (6.28) 



This expression is valid for (1 + z) eoeo ^ 1- 

6.4 Sparticle Effects - An Estimate 

We can expect effects from sparticles to show up in the PBH emission above T ~ 100 GeV (if 
we have a fairly low sparticle mass spectrum). To illustrate the behaviour of the diffuse flux 
from PBHs further, especially how SUSY scenarios might affect it, wc turn to studying the 
diffuse flux integral at peak photon energies. We begin by considering the low-temperature 
part of the integral (6.7) 

r*' 



^ oc ^J^Jl+..)%iTit.),E)dt.. 

« lO^s g^-i f (1 _^ zi)3/2/^(Ti, E)dzi < 



< 



^"''s/t''"' [^^ + - + ^'^'^1 ^^^"'^^ 

1025m-2sr-i 



5/2 



(1 H" -^max)^^ f-f{Ti, Ep^nl^) , (6.29) 



where t' is the time when the black hole temperature becomes some order of magnitude 

higher than the initial temperature T^, z^ax is an approximation for the highest contributing 
redshift and Ep),a,k ~ 100 MeV is the peak emission photon energy. This approach was 
chosen since the black hole spends most of its time near Mj, and during this time the flux 
does not change very much. It was assumed that Otot = = 1 for simplicity. In the 
last step we assumed Zmax S> z' which is valid for the range of PBHs we consider. The 
optical depth was also assumed to be significantly higher than 1 for z > z^ax and not to 
change significantly for lower redshifts. This should be sufficient for an order of magnitude 
estimate (within a factor 10 or so) in the observable photon energy range where one can 
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sec the furthest back in time. It will become clear later that only an upper limit such as 
this will be sufficient for our purposes. 

To get an estimate for the high temperature flux, the integral (6.7) is rewritten in terms 
of the black hole temperature T. First, we need the expression 



STTkG J 2$ 



°{Tr^-T-^) + {l + zr'^' 



-2/3 



(6.30) 



This expression neglects the dependence in $ on t, i.e. we choose to evaluate $ at some 
time, typically the initial time ti. The function $ generally varies less than a factor 10 over 
the PBH lifetime (at least for the PBHs we consider) , so our assumption will underestimate 
$ with at most a factor 10. When T ^ Ti and Zi ~ 10^^ the expression is simplified to 



Using this expression for 1 + we can write the high-temperature part of the integral (6.7) 
as 



""•^ oc 5x lO^-V-^sr-i ' ' 



dEo V 1 GeV 

/■To 

X 



X 



/ ° /(T)5/=^e-'"(^'"^i(^»^(T, E)T-^dT , (6.32) 
Jt' 



where T' = T{t') with t' as earlier and Tq is the temperature at the present time (or rather, 
the highest temperature the PBH has had by today in the case that the PBH has completely 
evaporated by today). Neglecting optical depth, and realising that the redshift at photon 
emission remains fairly constant, see Eq. (6.31), we may simplify the integral to 



dJ 
dEo 



CX 5 X 1033m-2sr-l (j^) M^ref , Epeak)T-ff X 

Jt' 

xl033m-sr-(3A_j ^/(T,)^/3 



3-/3 



ref 

X f-riTref , Epeak) ■ (6.33) 

It was here assumed that the instantaneous flux oc . This is supported by simulation 
data, with 1 < /? < 2 (see e.g. Figure 8.1 or Ref. [4]). We thus evaluate at some reference 
temperature Tref and photon energy Epeak ~ 100 MeV (since we are interested in the peak 
flux). We can choose Tref = T' to simplify the expression somewhat. The assumption of 
constant redshift of emission is further supported by the fact that around the instantaneous 
peak energy, the flux does not vary signiflcantly with photon energy. Hence, the flux is not 
very sensitive to variations in z at these photon energies. The last step assumes Tq ^ T' 
which will be true in the cases we consider. 

We can now use these assumptions for the cases we consider. In these cases, the low- 
temperature contribution (6.29) will in fact be negligible compared to the contribution from 
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higher temperatures, so we need only consider the case T ^ Ti. For a black hole with mass 
Mj « 5 X 10^'* g we have 

Ti w 0.02 GeV , 
T' « 0.1 GeV, 

Zmax ~ 10^ (6.34) 

and inserting these in Eq. (6.33) with T^e/ = T' we obtain the requirement 

_^0.1^100o_i-3 < fWO^Tp, 100-3 ^ fWO^Tp, > IqS/O.I^IOO ^ (g_35) 

for the fliix from T > 100 GeV to be non-negligible. Here f^^^'^^ symbolically denotes "the 
flux between T = Ti and T = T2" . In the case of a black hole with M, w 3 x 10-^^ g we use 

Ti w 0.35 GeV, 
T' « 1 GeV , 

Zmax ~ 2 X 10^ , (6.36) 

and in the same manner as above, we obtain the requirement 

^^1^100^-3 < /loo^T,, 100-3 ^ /^loo^Tp, > io«/i^i°" . (6.37) 

From simulation data (see Figure 8.1 or Ref. [4]) wo know that the fluxes are approximately 

_^o.i^ioo ^ io23GeV-is-i, (6.38) 

« 1025GeV-is^i. (6.39) 

Inserting this above we find that the instantaneous flux from temperatures above 100 GeV 
should fulfiU 

fWO^Tp, > ^q32 Mi « 5 X 10^^ g , (6.40) 

fWO^Tpi > 1031 GeV-^s-i, Mi w 3 X lO^^ g , (6.41) 

for the observed flux at peak energies to be significantly affected. These approximations 
should be valid to within about one or two factors of ten, taking into account also optical 
depth effects. The values agree well with MacGibbon's estimate that the increase in degrees 
of freedoms for T > 100 GeV would have to be of the order 10^ to significantly affect the 
lifetime emission of a PBH [5] . 

In the preceding, we have looked at the dominant peak emission, with emitted photon 
energy ^ 100 MeV. There is of course the possibility that the diffuse fiux coming from higher 
emitted energies may be affected. This will be explored further in the Results Section. 



6.5 Applied Model 
6.5.1 Initial Mass Function 

Green and Liddlc have shown [28] that a reasonable first approximation is to use a PBH 
distribution created at t ~ with a single initial mass M° and volume number density rij. 
This gives the initial mass function 

^ = n,<5(M, - MO) = no{l + z,)H{M, - Mf) , (6.42) 

where no is the present number density and Zi is the redshift of PBH formation. An 
important note to make is that this definition of the present number density is just the 
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diluted initial density, hence the effect of evaporation is not considered (i.e. no PBHs 
disappear). The value of no can be constrained using the diffuse gamma-ray flux from PBHs. 
It can also be independently constrained by requiring ClpBH < 1- Using the formulae of 
Section 2.1.4 and Ref. [18], one finds that the value at the epoch of PBH formation is 

where pc = 3H^ / (SttG) is the critical energy density. The value of Op^j^ is also an upper 
limit on flpBH at later times. Apparently, the value of no can be significantly greater than 
1 pc~^ for PBHs emitting particles in non-negligible amounts (i.e. with sufficiently small 
initial masses) without the Universe being overdosed. 
We will examine two "limiting" cases for M?, 

• a PBH expiring today 

• a PBH expiring when (the thinnest part of) the optical depth dropped below unity 
{z ~ 900) 

These limits correspond to the range of PBHs contributing to the DGB today, which have 
initial masses approximately in the range 3 x lO^'' - 5 x 10^"' g, corresponding to initial 
temperatures approximately in the range 0.02 — 0.35 GcV. We will thus use M° « 5 x lO^** g 
and M° w 3 x lO^'' g respectively for the two cases listed above. 

We will also look at some "intermediate" mass PBHs in the range 4 x 10^"^ g < Mi < 
1 X 10^** g, since the flux from these might possibly show a difference in the high-energy part 
of the diffuse photon flux. 

6.5.2 Cosmological Pcirameters 

We will use cosmological parameters according to WMAP results [67], except for flroT = 
1.0, Zrec = 1300 (this deviation from WMAP results is still within uncertainties, and does 
not affect our conclusions). In summary, we use 



^0 


= 2.725 K 


ho 


= 0.71 


Qtot 


= 1.0 


CIm 


= 0.28 




= 0.72 


fib 


= 0.044 




= 20 


^rec 


= 1300 



(6.44) 



and thus we have the Zdziarski-Svensson parameters 



T2.7 


« 1.01 


hdo 


= 1.42 


flo.i 


= 0.44 


Oo 


w 4.60 X 10"^° 



(6.45) 
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Distribution ot optical depth 




Figure 6.1. Unit optical depth boundciry using results in [71]. Cosmological parameters are 
Qtot = 1.0, fit = 0.044, ft = 0.71, Zri = 20, Zrec = 1300, To = 2.725 K according to WMAP 
results [67]. 



6.5.3 Optical Depth 

The optical depth was evaluated using the formulae given in Section 6.3. With the pa- 
rameters in the preceding Section, the optical depth then takes the form shown in Figure 
6.1. 



6.6 GLAST 

The Gamma-ray Large Area Space Telescope (GLAST) due for launch in 2006 will improve 
measurements of the diffuse gamma-ray background. GLAST contains two instruments, the 
Large Area Telescope (LAT) observing the gamma-ray background and the GLAST Burst 
Monitor (GBIVI) searching for gamma-ray bursts. 

A summary of the specifications of the LAT is shown in Tabic 6.1 along with comparisons 
to EGRET. The GBM specifications are reported in Table 6.2 with comparisons to BATSE. 
For a fuller description of these characteristics, see Ref. [3] . 

It is clear that GLAST measurements will significantly improve our knowledge of the 
diffuse gamma-ray background as well as gamma-ray bursts and hence related phenomena 
such as PBH abundance and primordial density perturbation constraints. 
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Quantity 


LAT (minimum spec.) 


EGRET 


Enca'gy range 


20 M(^V - 300 GcV 


20 McV - 30 GcV 


Peak (^libctiNC area"' 


> 8000 cui- 


1500 cm- 


Field of view 


> 2 sr 


0.5 sr 


Angular resolution'^ 


< 3.5° (100 MeV) 

< 0.15° (>10 GeV) 


5.8° (100 MeV) 


Energy resolntion*^ 


< 10% 


10% 


Deadtime per event 


< 100 lis 


100 ms 


Source location determination'^ 


< 0.5' 


15' 


Point source sensitivity'' 


< 6 X lO-^cm-^s"^ 


~ IQ-'^cm-^s"^ 



After background rejection 

Single photon, 68% containment, on-axis 

1-fT, on-axis 

1-a radius, flux lO-^cm-^s"! (>100 McV), high |&| 
° > 100 MeV, at high for exposure of one- year all sky survey, photon spectral 
index -2 



Table 6.1. LAT specifications compared to EGRET. Reproduced from [3]. 



Quantity 


GBM (minimum spec.) 


BATSE 


En(^rg\- raiigc^ 


< 10 k(A' - > 2.") McV 


23 k>V - 10 McV 


Field of view 


all sky not occulted by the Earth 


4p sr 


Energy resolution** 


< 10% 


< 10% 


Deadtime per event 


< 15 /is 




Burst sensitivity'' 


< 0.5cm~^s~^ 


0.2cm-^s-i 


Alert GRB location'^ 


~ 15° 


~ 25° 


Final GRB location^ 


~3° 


1.7° 



^ 1-a, 0.1-1 MeV 

^ 50 - 300 kcV 

Calculated on-board; 1 second burst of 10 photons cm~^ s~^, 50 - 300 keV 
Final ground computed locations; 1 second burst of 10 photons cm~^ s~^ 
50 - 300 keV 



Table 6.2. GBM specifications compared to BATSE. Reproduced from [3] 



Chapter 7 

Simulation 



7.1 Pythia 

The black hole emission was simulated using the particle physics event generator Pythia 
6.212 [73]. Pythia is a frequently used event generator using the Lund model of string 
fragmentation for QCD particles. It produces results well within experimental constraints 
for known physics and includes various theoretical models for hypothetical physics. Specifi- 
cally, it can simulate sparticle processes for various SUSY models. Pythia assumes R-parity 
conservation. Only the lowest order sparticle pair production is included. SUSY decays of 
the top quark are included, but all other Standard Model particle decays are unaltered. 
Pythia takes MSSM/mSUGRA parameters and calculates masses, cross sections etc. us- 
ing analytical formulae. 

Since Pythia is designed for beam collisions, and not direct decay of particles, some 
modifications had to be made, according to instructions from the author Torbjorn Sjostrand. 
The changes made were in the function PYRESD. For details, see Appendix B. 

7.2 Convolution Methods 

We wish to perform the convolution of the direct emission spectra with the fragmentation 
functions. To do this, one could choose (at least) two different approaches: 

1. Evaluate the fragmentation functions fully for some initial-particle energies, fit a 
scaling-relation for all energies and then integrate (most likely numerically). 

2. Generate a direct emission spectrum of particle energies and let them fragment and 
decay to find the final spectra. 

The benefit of the first option is that we obtain an expression for all temperatures T of the 
black hole. One of the downsides is that there are as of now no good parameterisations of 
the fragmentation functions for general coloured particles [74]. One could probably find 
fits, but this would require considerable computer time. 

The benefit of the second option is that we do not need any knowledge of how the 
fragmentation functions behave. It is also a more direct representation of the actual physical 
process. We save computer time from this, but the downside is that we only get the 
instantaneous flux for a specific black hole temperature T. We need to perform a new 
simulation for each black hole temperature we are interested in. However, between mass 
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thresholds, the instantaneous flux is not affected by any new emitted degrees of freedom, 
and so should behave fairly predictably. 

For our purposes, the second option was chosen. This is also the choice made by 
MacGibbon-Webber [4, 5] and Barrau et al. [37]. 



7.3 Rejection Method 

To perform the simulation we need to generate particle energies according to the direct 
(Hawking) emission spectrum. To achieve this, we use the rejection method initially intro- 
duced by von Neumann in 1951. A variant of it consists of the following algorithm: 

Definition 1 (Rejection method) We wish to generate random numbers according to 
the probability density p{x) in the interval [a,b]. Let Pm > maxa<x<b p{x) ■ Let Z be the 
desired random variable. 

1. Generate a random variable X e U{a,b) 

2. Generate a random variable Y e U{0,Pm) 

3. IfY >p{X), start over. 

4- Set Z=X, which is the generated random number (i.e. when Y < p{X)). 

Theorem 1 The rejection method produces the random variable Z distributed according 
to p{x) in [a, b] . 

Proof: 



P{Z <r)= P{X < r\Y < p{X)) 



To use the rejection method, we need to find the probability density for the directly emitted 
spectrum. Let Qj be a random variable equal to the total energy of an emitted particle of 
species j with rest mass Hj. We can then write (assuming q> iJ.j) 



P{Qj <q)= 7 , (7.1) 



since this corresponds to the fraction of emitted particles with Qj < q. We use the direct 
emission rate, Eq.(3.4). This is the probability distribution, so the probability density is 

dNj / \ 

which one might have seen straight away. 



7.4 Simulation Algorithm 

To perform our simulation, we can now use the following algorithm for a black hole with 
temperature T: 



Simulation Algorithm 
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Randomly select with probability 1/n a particle species j from the emitted species 
{pi}i K times (this is to do the same amount of simulation on each species). 

Find a particle energy Qj using the rejection method with probability density fq.. 
The limits ^ij/T < Q^/T < 10 can be suitable. 

Create the particle in Pythia, let it fragment and decay the products. 

Keep track of the different final products and store them for each initial particle 
species separately. 

When finished generating particles, normalise the fluxes to the total flux, i.e. multiply 

the stored fluxes coming from species j by ajNj/hj^ where fij « K/n is the number 
of j particles created during the program run. 
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8.1 Instantaneous Photon Flux 

The instantaneous photon fliox was simulated for the Standard Model and mSUGRA Bench- 
mark B (see Appendix A) according to the algorithm in Section 7.4. We limit the analysis 
to this mSUGRA benchmark model for reasons explained further in Section 8.3.2. The 
integrated fluxes Nj were found by numerically integrating the direct emission spectra. 
Numerical Recipes routines [75] were used for numerical integration, random number gen- 
eration and splines fitting. The results from the simulations are shown in Figure 8.1 - 
8.5. 

The dominant peak at « 100 MeV comes from QCD fragmentation prochicts, sig- 
nificantly pion decay (7r° 27). The low-energy, log-linear part comes from electroweak 
particles, mainly bremsstrahlung. The high-energy kink is the direct photon emission. 

Comparing our simulation results with those of MacGibbon and Webber [4], we find 
that for temperatures in the range T 10—100 GeV our results for the flux around the 
peak region is of the order 20% less, the difference decreasing with increasing energy. This 
difference seems to come from the difference in QCD masses used (for instance, [4] uses 
mt = 50 GeV, « 0.6 GeV whereas we use = 175 GeV, mg = OGeV). There seems 
also to be some difference in the treatment of the QCD threshold and direct pion emission, 
giving a difference in flux around T = 0.3 GeV of the order 50%. This is within the limits 
of the uncertainty in particle model behaviour around Aqcd [4]. 

As can be seen in Figure 8.6, the instantaneous flux from a PBH in mSUGRA Bench- 
mark B is roughly a factor 4 higher than in the Standard Model at peak energies, when all 
particles are emitted with relativistic energies. Comparing this with the expected needed 
increase in flux at high energies to produce a signiflcant effect on the DGB found in Section 
6.4, we see that it is too small by several orders of magnitude. We thus expect that no 
observable effect should be present in the dominant part of the photon spectrum. This is 
explored further in the calculation of the diffuse flux. 

8.2 Integrated Photon Flux 

Since each Standard Model degree of freedom has a supersymmetric partner degree of free- 
dom, one might expect that the total flux from a PBH in a supersymmetric scenario would 
be roughly twice that of a Standard Model scenario (at least in the limit were all particles 
are emitted at relativistic energies). However, one also has to consider the difference in 
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Total photon flux 
Standard Model 




1 2 

log,„(E [GeV]) 



Figure 8.1. Total photon flux from black holes in the Standard Model. Jagged curves are 
simulation data and smooth curves the fitting. The table to the right gives the black hole 

temperature T. 
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Figure 8.2. Total photon flux from a T=100 GeV black hole in the Standard Model and 
mSUGRA Benchmark B. 
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Figure 8.3. Total photon flux from a T=200 GeV black hole in the Standard Model and 
mSUGRA Benchmark B. 
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Figure 8.4. Total photon flux from a T=500 GeV black hole in the Standard Model and 
mSUGRA Benchmark B. 
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Figure 8.5. Total photon flux from a T=2000 GeV black hole in the Standard Model and 
mSUGRA Benchmark B. 
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Figure 8.6. The relative flux of a PBH in mSUGRA Benchmark B compared to the Standard 
Model. Note that the fluctuations at high temperature lie within uncertainties. 
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direct emission flux for different particle spins. A rough approximation in the hmit that all 
emitted particles are relativistic could then be 

?'^^^^ » 1 + 3.8-3.9, (8.1) 

where aj denotes degrees of freedom of particle/sparticle species j, Ng. is the total direct 
flux, Eq. (3.4), of sparticle species i with spin Sj (similarly the denominator is the Standard 
Model particle flux). Wc sum of course only over emitted particles that might eventually 
produce photons. Comparing this somewhat naive prediction in the relativistic limit, with 
our results for mSUGRA Benchmark B, T=4000 GeV: 

N^,susY ~ 1.97X lO^i s^i, 
N^,SM ~ 5.65xl03°s-S (8.2) 

one finds 

,susY ~ 3.5iV^,sM , (8.3) 
which given the roughness of the approximation is in reasonably good agreement. 

8.3 Photon Spectra from PBHs 

8.3.1 Point Source Flux 

The likelihood that one should happen to locate and measure the photon flux from a single 
expiring PBH seems low. On the other hand, if we were to observe the final emission stages 
of a PBH, we might easily impose limits on allowed particle models since here we have a 
significant difference in flux for the peak energy, as can be seen in Figure 8.7. Note that 
the time scale is arbitrary. There is a clear quantitative and qualitative difference between 
the Standard Model and mSUGRA Benchmark B curves. Measuring this curve would allow 
conclusions regarding the mass spectrum of the true particle model. 

8.3.2 Diffuse Flux 

We limit ourselves in the analysis to the mSUGRA Benchmark Model B. This is because 
Model B has the lowest mass spectrum of the benchmark models considered, hence poten- 
tially the greatest effect on the photon spectra. For details on the mSUGRA Benchmark 
Models, sec Appendix A and Ref. [8]. Since one does not expect the effect of sparticlcs to 
be highly significant, an upper limit such as this should suffice for the diffuse gamma-ray 
flux. 

To find the diffuse photon flux, the numerical data for the instantaneous photon flux 
obtained using Pythia was fitted and the integral (6.7) numerically integrated. The initial 
mass function and cosmological parameters used are described in Section 6.5. 

The fitting of the instantaneous flux was done in a number of steps. The functional form 
chosen was 

lGeV-is-1 ; fiin{E,T) X {E) + 
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Figure 8.7. Point source photon flux during the final days of PBH evaporation for the 
Standard Model and mSUGRA Benchmark B. Note that the time scale is ajrbitraxy. 

+Uoiy{E,T) X w':^:°,^^\e) , (8.4) 
where fun is a linear part, fpoiy a higher order polynomial and 



log 



10 



LGeVy 



-b 



(8.5) 



6 being the Heaviside step function. The upper limit \ogiQ{E/lGeV) = —2.25 is the 
approximate point where the polynomial part starts to dominate over the linear part. Emax 
is the approximate energy where the flux has dropped significantly. This will coincide with 
the point where the direct photon emission dominates. The polynomials have the form 



/H„(i;,T)«^p„(T) 

6 

UolyiE,T)^J2cr,iT) 



E 



E 



For T > 2000 the coefficients exhibit a limiting behaviour, approximately parameterisable 



as 



Cn ~ d„i In 



(8.8) 
(8.9) 



8.3. Photon Spectra from PBHs 



47 



n 




4^ 





0.00052 


26.91986 


1 


0.00015 


-0.38854 



Table 8.1. Coefficients q^^ for the Standard Model. 



n 




"nO 





0.47408 


24.54293 


1 


0.03341 


-0.94904 


2 


0.02473 


-0.40625 


3 


0.01237 


-0.09398 


4 


-0.00034 


-0.00562 


5 


-0.00127 


0.01257 


6 


0.00026 


-0.00241 



Table 8.2. Coefficients dfif for the Standard Model. 



n 


Inl 







0.00047 


27.35710 


1 


0.00012 


-0.46435 



Table 8.3. Coefficients q^^ for mSUGRA Benchmark B. 



n 


"nl 


d^ 





0.49768 


24.95461 


1 


0.05629 


- 1.23295 


2 


0.01511 


- 0.36052 


3 


0.00644 


- 0.00381 


4 


0.00409 


- 0.04291 


5 


-0.00040 


- 0.00150 


6 


-0.00016 


- 0.00225 



Table 8.4. Coefficients d?^ for mSUGRA Benchmark B. 
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Figure 8.8. Diffuse photon flux from a delta-function PBH mass spectrum (see Section 6.5), 
with Mi = 5.125 x lO" g. This corresponds to Thh Ri 4 x 10* GeV, zf-^ ss 5 x lO^'' and 
2? Ri6xl0^'' (see Eq. (2.11)-(2.12)). The normalisation is arbitrary, i.e. the flux is normalised 
to the present-day PBH number density no as defined in Section 6.5. The COMPTEL [62] 
and EGRET [61] values are of course absolute flux values, not normalised to PBH density. 
The apparent discontinuities before and close to the drop at Eq ~ 1 GeV are artefacts from 
the numerical treatment and should not be considered physical. 



When numerically integrating the diffuse photon flux, the coefficients dn, Qn were de- 
termined using cubic splines on data from least squares fits on the numerical data for the 
instantaneous flux. For temperatures T > 4000 GeV for which no numerical data existed, 
the Eqs. (8.8)-(8.9) were used. This approach proved reliable up to T w 5000 GeV. For 
higher temperatures, the temperature was approximated to T = 5000 GeV. This produces 
a negligible error in the diffuse flux. 

The value of Emax, the highest energy that gives a significant fiux, was determined using 
cubic splines on simulation data. Above Emax, the direct photon flux was still included. 

Fitting of the instantaneous flux using Gauss-Seidel to various non-linear functions with 
a Boltzmann factor, a functional form one might expect would be successful, were also 
attempted but did not yield good accuracy. 

The results from the numerical integration is shown in Figs. 8.8 - 8.18. Note that the 
Figures have an arbitrary normalisation, given by the present-day PBH number density no 
as defined in Section 6.5. The peak in the diffuse spectrum is located at roughly i?o,peafe ~ 
500 keV. This peak comes from the dominant photon production due to QCD fragmentation 
(primarily pion decay), with energy at emission Ep^ak ~ 100 MeV. The high-energy part 
comes mainly from direct photon emission. 

As can be seen in Fig. 8.8, the effect from sparticles in the dominant part of the 
spectrum is completely negligible. For Mi « 5 x 10"'^^ g, there is a ~ 1 — 5% higher flux for 
mSUGRA Benchmark B for high energy photons (Figure 8.9), completely within numerical 
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Diffuse photon flux, M. « 5 x 1 g 
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Figure 8.9. Diffuse photon flux from a delta-function PBH mass spectrum (see Section 6.5), 
with Mi = 5.125 X lO"*^^ g. The normalisation is arbitrary, i.e. the flux is normalised to the 
present-day PBH number density no as defined in Section 6.5. Zoom of high-energy part of 
spectrum. The mSUGRA Benchmark B flux is of the order 1 — 5% higher than that for the 
Standard Model. 
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Figure 8.10. Diffuse photon flux from a delta-function PBH mass spectrum (see Section 
6.5), with Mi 3 X 10^-^ g. This corresponds to Tbh « 2 x lO'^ GcV, zf^^ 3 X 10^* and 
ai4xlO'^* (see Eq. (2.11)-(2.12)). The normalisation is arbitrary, i.e. the flux is normalised 
to the present-day PBH rmniber density no as defined in Section 6.5. The COMPTEL [62] 
and EGRET [61] values are of course absolute flux values, not normalised to PBH density. 
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Figure 8.11. Diffuse photon flux from a delta-function PBH mass spectrum (see Section 
6.5), with Mi PS 4.5 x lO^^ g. This corresponds to T^h 1 x 10^ GeV, zf^ « 2 x 10^8 
and zfraZY- IQp^ (see Eq. (2.11)-(2.12)). The normahsation is arbitrary, i.e. the flux is 

normalised to the present-day PBH number density no as defined in Section 6.5. The EGRET 
[61] values are of course absolute flux values, not normalised to PBH density. 



Chapter 8. Results 



Diffuse photon flux, M. « 4.5 x 1 g 



Standard Model 
mSUGRA Benchmark B 




9.1 9.15 
log,„(E„)[eV] 



Figure 8.12. Diffuse photon flux from a delta-function PBH mass spectrum (see Section 
6.5), with Mi Si 4.5 X 10^^ g. The normalisation is arbitrary, i.e. the flux is normalised to 
the present-day PBH number density no as defined in Section 6.5. Zoom of high-energy part 
of spectrum. The Standard Model flux is of the order 10% higher than that for mSUGRA 
Benchmark B. 
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Figure 8.13. Diffuse photon flux from a delta-function PBH mass spectrum (see Section 
6.5), with Afj ^ 5.5 X lO^^ g. This corresponds to Trh » 1 x lO'^ GeV, 2 X 10^* (see Eq. 
(2.11)-(2.12)). The normaUsation is arbitrary, i.e. the fiux is normalised to the present-day 
PBH number density no as defined in Section 6.5. The EGRET [61] values are of course 
absolute flux values, not normalised to PBH density. 
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Figure 8.14. Diffuse plioton flux from a delta-function PBH mass spectrum (see Section 
6.5), witli Mi 5.5 X 10^^ g. The normalisation is arbitrary, i.e. the flux is normalised to 
the present-day PBH number density no as defined in Section 6.5. Zoom of high-energy part 
of spectrum. The Standard Model flux is of the order 20% higher than that for mSUGRA 
Benchmark B. 
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Figure 8.15. Diffuse piioton flux from a delta-function PBH mass spectrum (sec Section 6.5), 
with Mi Si 6 X 10^^ g. This corresponds to Trh ~ 1 x f GcV and as 2 X 10^* (see Eq. 
(2.11)-(2.12)). The normaUsation is arbitrary, i.e. the flux is normalised to the present-day 
PBH rmmber density no as defined in Section 6.5. The EGRET [61] values are of course 
absolute flux values, not normalised to PBH density. 
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Figure 8.16. Diffuse photon flux from a delta-function PBH mass spectrum (see Section 
6.5), with Mi 6 X 10^'^ g. The normalisation is arbitrary, i.e. the flux is normalised to the 
present-day PBH number density no as defined in Section 6.5. Zoom of high-energy part of 
spectrum. The Standard Model flux is of the order 25 — 30% higher than that for mSUGRA 
Benchmark B. 
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Figure 8.17. Diffuse plioton flux from a delta-function PBH mass spectrum (see Section 6.5), 
with Mi ^7 X 10^3 g. This corresponds to Trh « 1 x 10^ GeV and fs 2 x 10^* (see Eq. 
(2.11)-(2.12)). The normalisation is arbitrary, i.e. the flux is normalised to the present-day 
PBH number density no as defined in Section 6.5. The Standard Model flux is of the order 
40 - 50% higher than that for mSUGRA Benchmark B. 



Chapter 8. Results 




Figure 8.18. Diffuse photon flux from a delta-function PBH mass spectrum (see Section 6.5), 
with 8 X 10" g. This corresponds to Trh ~ 1 X 10^ GeV and 2^ 2 x 10^* (see Eq. 

(2.11)-(2.12)). The normalisation is arbitrary, i.e. the flux is normalised to the present-day 
PBH number density no as defined in Section 6.5. The Standard Model flux is of the order 
25 - 30% higher than that for mSUGRA Benchmark B. 



8.4- Other Particles 
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uncertainties. Comparing the diffuse flux curve to that of MacGibbon and Carr [59], one 
finds that we do not have a bump at Eq 100 MeV, but rather a drop at Eq 500 — 
1000 MeV. This difference seems to depend on the difference in initial mass function (they 
assume a Carr initial mass function) and cosmological parameters. The bump, or drop, 
comes from the optical depth. The flux at Eq > 1 GeV comes from the high-temperature 
stage of emission with z 0. Note that the apparent discontinuities close to the drop are 
artefacts from the numerical treatment, and should not be taken as physical. 

Figure 8.10 shows that also for Mj « 3 x 10^^ g the supersymmetric contributions are 
negligible. Furthermore, the optical depth effectively cuts of all photons with Eq > 60 MeV. 

Turning to other initial masses (Figs. 8.11-8.18), one sees that in the high-energy part 
of the spectrum, a bit counter-intuitive perhaps, the Standard Model flux is greater. In 
this part of the spectrum, it turns out that for a range of initial masses (4 x 10^^ g < < 
1 X lO^"* g) the faster increase in black hole temperature associated with a SUSY particle 
spectrum is the dominating effect (over difference in instantaneous flux). This is due to 
the fact that the SUSY black hole spends less time "at each temperature" . The lower the 
sparticle mass spectrum is, the broader is the photon energy range in which this effect 
will be present. For a higher mass spectrum than in mSUGRA Benchmark B, one would 
then expect a similar effect, but "beginning" at higher photon energy. The magnitude of 
the difference will be dominated by the optical depth, since the effect comes from the late 
stages of PBH emission (previous emission needs to be attenuated not to dominate the late 
stage emission). Because of this, PBHs that reach the final stages of emission close to when 
the universe becomes thin to the photon energy range in question will show the greatest 
difference, since previous emission is highly attenuated. The results show that the effect 
occurs for Eo > 1 GeV, and the greatest difference is for Mj w 7 x 10^^ g where the Standard 
Model flux is ^ 40 — 50% greater than that for mSUGRA Benchmark B. For smaller initial 
masses, more of the interesting photon energy range is cut off by the optical depth. For 
larger initial masses the effect is still present, but quickly becomes dominated by previous 
emission, turning over to a Standard Model flux ^ 1 — 5% greater than the mSUGRA 
Benchmark B flux for Mj 5 x lO^'^g. The total effect on the observed spectrum will of 
course be dependent on the shape of the PBH initial mass function, but should in general 
be much less than 50%. It should be pointed out that the treatment of the function $ in 
the mass loss rate was not very exact close to mass thresholds (step function), which might 
affect the results slightly, although not significantly. 

Additionally, normalising the spectrum in Figure 8.8 to that measured by COMPTEL 
[62] and EGRET [61], one sees that the diffuse PBH photon flux at energies Eq > 1 GeV 
is some orders of magnitude less than the measured one. Hence, we do not expect to 
discriminate between the Standard Model and MSSM particle models by measLuing the 
diffuse photon flux from PBHs, at least if PBHs with initial mass Mi 5 x 10^^ g dominate. 

8.4 Other Particles 

The instantaneous flux of other astrophysically stable particles show the same general 
behaviour as the photon flux as regards differences between the Standard Model and 
mSUGRA. This was explicitly checked from the simulation data obtained. We will not 
consider these particles further in this work. 
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Chapter 9 

Conclusions 



Emission from PBHs remains an interesting testing ground for high-energy physics and the 
conditions in the early Universe. In these considerations, the effects of supersymmetric 
particle models on the gamma-ray flux from PBHs were studied. 

The numerical data for the instantaneous photon flux from PBHs found agrees quite 
well with previous works, although some differences exist due to some difference in model 
assumptions. The instantaneous flux in mSUGRA models was found to be roughly a factor 
4 greater than in the Standard Model, for peak photon energies in the regime were all 
emitted particles are relativistic. 

The instantaneous flux of other astrophysically stable particles was found to exhibit the 
same general behaviour as the photon flux. 

For the point source flux, the shape of the flux at peak energies as a function of time 
will be quantitatively and qualitatively different for the Standard Model and mSUGRA, 
allowing limits to be imposed on the true particle model's mass spectrum. Albeit a far- 
fetched prospect to be able to detect the emission from a PBH in its final year to days, it 
would provide valuable information on a variety of physics, particularly the mass spectrum 
of the true particle model. 

The results show that no observable difference should be present in the dominant part of 
the diffuse gamma-ray flux. In the high-energy range Eq > 1 GeV, a PBH in the Standard 
Model produces a diffuse flux that is typically much less than ~ 50% greater than for a 
PBH in mSUGRA Benchmark B. The exact difference will be dependent on the shape of 
the PBH initial mass function. For other SUSY models, the photon energy range where 
this difference is present will start at higher or lower energies for heavier and lighter particle 
mass spectra respectively. In the context of DGB measurements, particularly the upcoming 
GLAST project, it seems unlikely although perhaps not impossible to see such an effect. A 
detailed study using a specific initial mass function is needed to settle the question definitely. 

The effect of possible photospheres would be to attenuate the low-energy (QED photo- 
sphere) and high-energy (QCD photosphere) flux. A detailed analysis is outside the scope 
of this work. 
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Appendix A 

mSUGRA Benchmark Models 



A.l B 

The parameters of mSUGRA Benchmark B according to Ref. [8] is shown in Table A.l. 
The resultant mass spectrum of non-Standard Model particles is shown in Table A. 2. 



A.2 A, C-M 

For parameters, see Ref. [8]. The mass spectra can easily be calculated using e.g. Pythia 
by giving the parameters as input. 



Parameter 


Value 


mi/2 


255 GcV 


Too 


102 GcV 


tan(/3) 


10 


sgn(M) 


+1 


A 






Table A.l. Parameters for mSUGRA Benchmark B 
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Appendix A. mSUGRA Benchmark Models 



Partick 




Alass \CicY] 


x¥ 




101.73 






111.64 


4 




148.49 




148.49 
148.89 


Xi 




1 on Afi 


A2 




191.24 


^IJ,,Lj I'll 


,R 


200.14 


l'e,L, t^e, 


R 


200.14 


4 


R 


205.58 




215.41 




215.41 






223.96 






366.75 


x\ 




386.69 






387.30 


AO 




408.54 






408.80 






416.11 






451.26 


6i,6i 




557.60 






581.95 


dR,_dL 




581.95 






582.97 


cr, cl 




583.11 


UR,Ul 




583.11 


Cl,Cr 




601.13 


UL,UR 




601.13 






605.79 


cIljJr 




605.79 


t2,t2 




631.22 


9 




681.31 



Table A.2. Sparticle/non-standard particle mass spectrum in mSUGRA Benchmark B 



Appendix B 

Pythia modifications 



The changes made in Pythia 6.212 were in the function PYRESD, and follow below. 
Line 13058: 

C . . . Order incoming partons and outgoing resonances . 
C***old IF ( JTMAX . EQ . 2 . AND . ISUB . NE . . AND . MSTP (47) . GE . 1 . AND . 

IF(JTMAX.EQ.2.AND.MSTP(47) .GE. LAND. ! new line 

Line 13190 and 13210: 

C... Store incoming and outgoing momenta, with random rotation to 
C... avoid accidental zeroes in HA expressions. 
C IF(ISUB.NE.O) THEN ! removed (made comment) 

DO 470 I=IMIN,IMAX 
( . . . etc . . . ) 
500 CONTINUE 
C ENDIF ! removed (made comment) 

Line 13233 and 13246: 

C. . .Calculate four-products. 

C IF(ISUB.NE.O) THEN ! removed (made comment) 

DO 540 1=1,2 
( . . . etc . . . ) 
560 CONTINUE 
C ENDIF ! removed (made comment) 

ENDIF 
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